Clustering of Fermi particles with arbitrary spin 
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A single Z-shell model is investigated for a system of fermions of spin s and an attractive s-wave, 
spin channel independent, interaction. The spectra and eigenvectors are determined exactly for 
different I, s values and particle numbers N. As a generalization of Cooper pairing it is shown that 
when N = /i(2s -\- 1), jj, — 1,2, ... ,21 + 1, the ground state consists of clusters of (2s + 1) particles. 
The relevance of the results for more general situations including the homogeneous system is briefly 
discussed. 
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PACS numbers; 32.80.Pj,05.30.Fk,74.20.Fg 

The experimental realizations of hi^ly degenerate 
atomic Fermi gases confined in traps [ll l2j S IJj 1^ Ifl t3| 
have increased the activity in the theoretical investiga- 
tion of such systems, which had already started earlier 
in view of the hope of future experiments. In particular 
many aspects of the possible transition to the superfluid 
phase have been studied in the case of the attractive in- 
teraction between the particles 0, Ip. Eol ITll ll^ ll^ Q, 

IIlIIlIl2bj^jj,li3,IIllil lilliMlalar For 

instance, the critical temperature for the Cooper pairing 
of atoms with spin values higher than 1/2 has been inves- 
tigated |il2i, J^J based on the Gorkov-Melik-Barkhudarov 
approach to calculate corrections to the BCS theory due 
to induced interaction |28| . The result is that the crit- 
ical temperature of the BCS theory is multiplied by 
(4g)i^/3-i^ where v denotes the number of fermion species 
(v = 2s+l) T?]. Accordingly, if s = 1/2 (the case treated 
by Gorkov-Melik-Barkhudarov) the critical temperature 
is decreasing, while it is increasing for s > 3/2. In the 
special case s = 1 no such correction exists |l2l |. 

The aim of the present paper is to propose and demon- 
strate that the superfluid ground state can be different 
from that given by the BCS pairing theory if the spin of 
the atoms s > 1/2. Namely, we suggest that instead of 
Cooper pairs the particles can create clusters of zero total 
spin value containing 2s + 1 particles. This can occur in 
optical traps where the spins of the particles can freely ro- 
tate (Alternatively, one can think about 2s -I- 1 hyperfine 
states. In the following we use always the terminology of 
spin). We do not attempt to make any comparison with 
experimental results in this paper and choose the pos- 
sible simplest model to make the presentation free from 
approximations. We hope that in future experiments the 
clustering effect can be revealed. We treat the case when 
the number of species (2s -I- 1) is even. It is straight- 
forward to extend our considerations for odd number of 
species. Note that in the spin-1 Fermi superfluid the pair 



correlation has been investigated in detail by Modawi and 
Leggett 

We assume that the trap potential is spherically sym- 
metric. The investigation will be done at zero tem- 
perature and is restricted to the subspace of a given 
Z-shell filled with particles of number (0 < iV < 
(2Z-|-l)(2s-f 1)). The interaction between the particles in 
the open shell is described by a spin independent contact 
potential representing s-wave scattering with a negative 
scattering length: 
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The one particle normalized wave functions in the open 
shell are given by 

«'«,i,m„s,m,(r,1?,(^,(T) =i?„,i(^)^4(^''^)Xm.M' (2) 

Here cr is a discrete spin variable, which can take 
(2s -f 1) different values. Spin eigenfunctions are ortho- 
normalized according to 



(3) 



The functions lO are eigenfunctions of the one-particle 
Hamiltonian, which contains besides the trap potential 
also the average field of the closed shells. The problem of 
diagonalizing the operator |^ can be solved exactly on 
the fixed basis 101. 

Let us denote by ami.m, the operator which annihi- 
lates a particle with quantum numbers (n,l,mi,s,ms). In 
second quantization ^ can be written as 



Hint=hEo/n= ^ 

7ni ,m2 ,7713 ,7714 ^1:^2 



m2 ,7773 ^'>TT'4 



mi ,1/1 ^m2 ,1/2 ^^114,1^2 "-ms ,1/1 , 
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where Eq is the characteristic energy 
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and the coefScient / is given in terms of the Wigner-3j 
symbols. 

The Hamihonian can be converted into 
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where the abbreviation [a] = (2a + 1) is used and the 
scalar "square" 
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of the irreducible tensor operators 

I s 



m=~l u=~s 



I I L 
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i/'^m — M,i' 
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has been introduced. N stands for the particle num- 
ber operator. We note that in cases of / = 1 and 
I — 2 the Hamiltonian can be expressed in terms 
of the operators N, C2(5'C/(3)), C2{SO{?,)) and iV, 
C2(5'C/(5)), C2(S'0(5)), respectively, where C2(G) de- 
notes the quadratic Casimir operator of the group G [i^ . 

We have calculated the full exact spectra and the eigen- 
vectors of h for s = 1/2, / = 1,2,3 and for s = 3/2, 
I — 1,2 using some Mathematica and Fortran pro- 
grams. Unfortunatelly, the size of the Hilbert-space (in 
the Lz — Q Sz ~ Q or Sz — 1/2 subspace too) grows 
drastically. For higher s we also have performed calcu- 
lations for Z = 1, s = 5/2,7/2 in the L = subspace 
for even N . The Young tableau technique has also been 
used to analyse the model. The details will be published 
elsewhere 2B\. Results for the ground state energies up 
to the half filled shell can be found in Tables and 
nil Above the half filling energies and eigenvectors can 
be obtained by particle-s-^hole transformation. If one ex- 
act eigenstate has N particle in the /, s-shell with energy 
E{N) then there is an other exact eigenstate for particle 
number N' = Nt~N {Nt = (2s + 1)(2/ + 1)) with energy 

E{Nt~N) = E{N) + {2l + l)s{Nt/4:-N/2), N < Nt/2 

It is enlightening to discuss the structure of the ground 
state for even number of particles in coordinate represen- 
tation as well. Let us start with the two particle states. 
The wave function is as follows: 



V'(l,2) = ^(ri,r2)^^+^X(^1,^2) 



(9) 



The function ip is symmetric and is the eigenfunction of 
the angular momentum operator with L = 0, while the 
spin functions are antisymmetric. The latter is ensured 
if the total spin S takes the values 5 = 0, 2, . . . , 2s — 1, 
which makes the ground state 2s times degenerate. The 
degeneracy is lifted and a spin multiplet arises, if the 
s-wave scattering lengths were different in different spin 
channels. In case of s = 1/2 it is found that the exact 
ground states are "pair states" for = 2, 4, . . . , 2(2/ -fl), 
i.e., their wave functions are the antisymmetrized prod- 
ucts of singlet pair wave functions 



*(1,2,...,7V) 



N/2 



(10) 



i=l 



The operator A stands for antisymmetrization. The sit- 
uation is different when s > 1/2. Namely, instead of pair 
states we find that 2s + 1 particle clusters are created 
when N = fi{2s + 1), fj. = 1, 2, . . . , 2/ -h 1. The wave 
function reads as 

^'(1,2,...,7V) (^(r,,r,+i,...,r,+2«) 

i 

^X(0-i,0-i+l, . . . ,0-i+2s) (11) 

The prime on the product sign in Eq. (|ll|l means that i 
extends to the values as follows: 



i = (2s + l)fc + l, 



A; = 0, 1, . . . , /i — 1 



(12) 



denotes the singlet spin function representing the 
Slater determinant of the (2s+l) one particle spin fun- 
tions Xm, 



^Xi<^l:(^2, ■ ■ ■ , CT2s+l) = 



1 



V(2s + 1)! 

X-s(cri) X-s(o-2) ••■ X-s{cr2s+l) 
X-s+l{cri) X-s+l{<J2) ■■■ X-s+l{cr2s+l) 



Xsi(Ti) Xs(cr2) 



Xs{<J2s+l) 



(13) 



In Eq. (|11|) (/3 is a completely symmetric function for the 
exchange of any pair of particles. Furthermore, it is an 
eigenfunction of the orbital angular momentum operator 
of the 2s + 1 particles with eigenvalue zero. 

Eq. Hll|l means that the states can be created by re- 
peated application of a creation operator to the vacuum 

Q2S+1 ) We have derived this operator, but we do 



not present here in general form. For s = 1/2 the opera- 



tor Q2 '' is the pair creation operator 
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FIG. 1: The energy levels of the dimensionless Hamilton 
operator h in case of Z = 2, s = 3/2. The dotted line (drawn 
everywhere to lead the eye) at even A'^ gives the average energy 
calculated using the pair wave function. The insertion shows 
the binding energy per particle as a function of the number 
of particles. 



and for the ground state energy the simple relationship 
is obtained for s = 1/2: 

9/ + 1 

Eo = —^N, iV(even) < 2(21 + 1) (15) 

Furthermore, for iV = 2 Eq. H15|) gives also the ground 
state energy for any s values. 

One can show that in the single Z-shell model the clus- 
ter wave function ip in Eq. I|ll|l has the structure 

s+l/2 

(/?(ri, . . . ,r2s+i) = S* Y[ 'P(i"2j-i,r2j), (16) 
i=i 

where S denotes the operator for symmetrization and 
(/j(ri,r2) is the function defined in Eq. It can be 

expressed in terms of the one particle wave functions: 



(/3(ri,r2) = i?n,i(ri)i?„, 2(7-2) 



E (-l)'"™il(^^l)>^-™(^^2). (17) 



In Fig. n besides the exact levels for I = 2, s = 3/2 
the expectation value of the Hamiltonian h, Eq. |0J, 
with wave function (|10|l is also shown. As can be seen 
it provides the exact ground state energy only for N = 
2, (2s+l)(2Z+l)-2, (2s+l)(2Z+l). In the other cases the 
energy obtained with the help of ((Tn|) lies considerably 
higher than the ground state energy. This is especially 
true when N = ^(2s + 1), and the wave function l(TT|) 
applies. The binding energy per particle has maximum 
at these particle numbers. 

For understanding the origin of the ground state wave 
function Hll|l one has to realize that the relevant group 



of the model is SU{21 + 1) (g> SU{2s + 1). With the usual 
decomposition one gets states characterized by quan- 
tum numbers specifying the irreducible representations 
of SU{21 + 1) D SO{2l + 1) D S'0(3) and similarly for 
SU{2s + 1), where the relevant quantum number is the 
total spin related to the irreducible representations of 
SU{2) C SU{2s + 1). The seniority is associated with 
50(2/ + 1), while the orbital angular momentum L with 
50(3). The validity of the cluster wave function ^ can 
be best understood by studying first the structure of its 
spin function. For fJ- — I the spin function is obviously 
invariant under the group SU{2s + 1). It means that 
the generators of SU (2s + 1) commute with the operator 
Q2I+1 introduced in the text above (|14|l . Since ^ is ob- 
tained by repeated application of Q2S+1 it follows that 
is also invariant under SU{2s + 1) (consequently 5 = 0). 
The corresponding Young tableau consists of fi columns 
of length (2s -|- 1), which is known to specify a one di- 
mensional irreducible representation of SU{2s + 1) (see 
e.g. "soj). The Pauli principle to be fulfilled by requires 
that the corresponding Young tableau of SU (2Z + 1) must 
consist of /i rows, each of length (2s + 1). It contains for 
each fi a one dimensional representation of 50(2/ + 1), 
which leads to i = 0. The state is unambigously deter- 
mined by choosing this representation of SO{2l + 1). We 
emphasize that the symmetry of 4* is higher than just 
the rotational invariance in the coordinate space and in 
the spin space. As a matter of fact the symmetry is the 
highest possible for particle numbers N — /i(2s -I- 1). 

In summary, our results are in harmony with the ex- 
pectation that the ground state of a system of Fermi 
particles with attractive interaction is the possible most 
symmetric one with respect to the exchange of the space 
coordinates of the particles. This feature manifests itself 
spectacularly, when the number of particles is the mul- 
tiple of the number of species 2s -I- 1. Though in this 
paper the single /-shell model has been treated one can 
presume the validity of the conclusions more generally. 

As a generalization one can take into account more 
than one shell in constructing the cluster wave function 
Lp (configuration interaction) . Note that the totally anti- 
symmetric spin functions ^x, Eq. (|13|l . would not change. 
Another strategy would be to treat first the homogeneous 
system and apply a local density approximation for the 
trapped gas. The problem of the ground state of the ho- 
mogeneous system is, of course, interesting in itself. In 
this case tp is translationally invariant. When s = 1/2 
the wave function 5*, Eq. Ijlll) . which coincides now with 
Eq. (|10l) is the well known expression of the BCS ground 
state projected onto the N particle state and written in 
coordinate representation (see e.g. HJ)- For the spin-1 
Fermi gas, introduced by Modawi and Leggett [l3], we 
predict a ground state containing three particles, clusters 
of zero spin and with a symmetric iy9(ri, r2, rs) function 
in the expression (|ll|l . 
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TABLE I: Ground state energies of the dimensionless Hamil- 
ton operator h for I ^ 1 and for = 0, 2, . . . , 3(2s + l)/2. 



N 


2 


4 6 


8 10 


s = 1/2 


-5/4 


-5/2 




s = 3/2 


-5/4 


-65/14 -185/28 


-75/7 -375/28 



TABLE II: Ground state energies of the dimensionless Hamil- 
ton operator h for I = 2 and for A?" = 0, 2, . . . , 5(2s + l)/2. 



As an example let us assume that four particles with 
spin 3/2 are put outside the Fermi sea of noninteracting 
particles as a generalization of the original Cooper prob- 
lem for two particles (see for its discussion |^]). One 
has to emphasize that even that part of the wave func- 
tion, which depends on the space coordinates does not 
have the "pair structure" (|16|) in general, it is valid only 
within the single ^-shell model. One can use (|16|) as an 
ansatz, however, in a homogeneous system and compare 
it with 11U|) for the ground state. Now (^(ri,r2) is speci- 
fied by the requirement that it should satisfy the equation 
for a bound Cooper pair (because of the S function na- 
ture of the two particle interaction one needs the usual 
regularization (see e.g. .11])). The result of a simple, but 
somewhat lengthy calculation leads to: 



^pairs -^Acluster 



5A 



______ +0(i,) 

(18) 

where Epairs and E^ciuster are the expectation values of 
the Hamiltonian with the wave functions of form I|1U|) 
and (|11I16|) respectively, using the interaction iQJ. V 
denotes the volume of the system. It has been used that 
(^(ri,r2) = <;5(|ri — r2|) can be taken real. The function 
C(r) is defined as C(ri — r2) = / d^rs (/?(|ri — r3\)ip{\r3 — 
r2|). As expected the difference is of 0{1/V). Since 
if and C are dominantly positive the result shows that 
the wave function (fTT|l is energetically favored due to the 
fact that Hll|) is "more symmetric" in space coordinates 
than ((Tnt . To get an energy gain of 0(1) a more accurate 
ground state wave function is needed, which can account 
for a possible four particle bound state. To find it goes 
beyond the scope of the present paper and is a part of 
our planned investigations. 
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